DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCH 
MONTEREY,  CALIFORNI 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

THE 

EFFECTIVENESS  OF  HEAT  EXCHANGERS 
WITH  ONE  SHELL  PASS  AND 
FIVE  TUBE  PASSES 

by 

Lawrence  E.  Hess 

September  1985 

Th 

ssis 

Advisor:                Allan 

D. 

Kraus 

Approved  for  public  release;  distribution  is  unlimited 


1222B21 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Enter, d) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.     REPORT  NUMBER 


2.  GOVT   ACCESSION  NO 


3.     RECIPIENT'S  CATALOG   NUMBER 


4.     TITLE  (and  Subtitle) 

The  Effectiveness  of  Heat  Exchangers 
with  One  Shell  Pass  and  Five   Tube 
Passes 


5.     TYPE  OF   REPORT  &   PERIOD  COVERED 

Master's  Thesis; 
September  1985 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  author^; 

Lawrence  E.  Hess 


8.     CONTRACT  OR  GRANT  NUMBERfsJ 


9.     PERFORMING  ORGANIZATION   NAME   ANO  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93943-5100 


10.     PROGRAM   ELEMENT.  PROJECT,   TASK 
AREA  4   WORK  UNIT  NUMBERS 


11.     CONTROLLING  OFFICE   NAME   AND   ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93943-5100 


12.   REPORT  DATE 

September  1985 


13.      NUMBER  OF   PAGES 


NUMI 

121 


U.     MONITORING   AGENCY   NAME   4    ADDRESSf//  different  from   Controlling  Office) 


15.     SECURITY   CLASS,   (ol  Ihlt   report) 

Unclassified 


15a.     DECLASSIFICATION     DOWNGRADING 

SCHEDULE  i 


16.     DISTRIBUTION    ST  A  T  EM  EN  T  (ol  this  Report) 

Approved  for  public  release;  distribution  is  unlimited 


17.     DISTRIBUTION   STATEMENT   (ol  the   abstract  entered  In  Block  20,    If  different   from  Report) 


18.     SUPPLEMENTARY   NOTES 


19.     KEY   WORDS  ( Continue  on  reverse  aide  It  necessary   and  Identity  by  block  number) 


Heat  Exchanger 
Effectiveness  N 
Five  Tube  Pass 


tu 


Method 


20.      ABSTRACT  (Continue  on  reverse   side   It  necessary  and  identify  by  block  number) 

Heat  exchangers  with  one  shell  pass  and  n  tube  passes  are 
often  referred  to  as  1-n  exchangers.   The  heat  transfer 
literature  contains  many  references  to  studies  of  1-n 
exchangers  when  n  is  even  but  apparently,  other  than  a  single 
study  and  some  work  pertaining  to  the  1-3  exchanger  at  the 
Naval  Postgraduate  School,  little  work  has  been  done  with 
respect  to  the  1-n  exchanger  when  n  is  odd.   This  thesis  expands 
upon  the  study  of  1-n  exchangers  with  n  being  odd. 


DD     ,    JAN    73     1473  EDITION   OF    1   NOV  65  IS  OBSOLETE 

S    N   0102-  LF-  0  14-  6601 


SECURITY   CLASSIFICATION   OF   THIS  PAGE  (When  Data  Entered) 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWhan  Dm*  Enffd) 


20.   (Continued) 

While  a  completely  closed  form  solution  was  found  to  be 
unfeasible,  a  polynomial  approximation  has  been  developed 
that  yields  the  effectiveness  (e)  of  the  two  possible  arrange- 
ments of  the  1-5  exchanger  as  a  function  of  the  capacity  rate 

ratio  (R)  and  the  number  of  transfer  units  (N.  ) .   This  will 

tu 
enable  the  analyst  to  consider  exchangers  where  the  inlet 

to  and  outlet  from  the  tubes  are  at  opposite  ends  of  the 

exchanger . 


S    N   0102-  LF-  014-  6601 


SECURITY  CLASSIFICATION  OF  THIS  PAGE(T»Ti»n  Dmtm  Enttrmd) 


Approved  for  public  release;  distribution  is  unlimited 


The  Effectiveness  of  Heat  Exchangers 
With  One  Shell  Pass  and 
Five  Tube  Passes 


by 


Lawrence  E.  Hess 
Lieutenant  Commander,  United  States  Navy 
B.S.,  The  Citadel,  1973 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
September  1985 


0.1 


ABSTRACT 


Heat  exchangers  with  one  shell  pass  and  n  tube  passes 
are  often  referred  to  as  1-n  exchangers.   The  heat  transfer 
literature  contains  many  references  to  studies  of  1-n 
exchangers  when  n  is  even  but  apparently,  other  than  a 
single  study  and  some  work  pertaining  to  the  1-3  exchanger 
at  the  Naval  Postgraduate  School,  little  work  has  been  done 
with  respect  to  the  1-n  exchanger  when  n  is  odd.   This 
thesis  expands  upon  the  study  of  1-n  exchangers  with  n  being 
odd. 

While  a  completely  closed  form  solution  was  found  to  be 
unfeasible,  a  polynomial  approximation  has  been  developed 
that  yields  the  effectiveness  (e)  of  the  two  possible 
arrangements  of  the  1-5  exchanger  as  a  function  of  the 
capacity  rate  ratio  (R)  and  the  number  of  transfer  units 
(N-t-u).   This  will  enable  the  analyst  to  consider  exchangers 
where  the  inlet  to  and  outlet  from  the  tubes  are  at  opposite 
ends  of  the  exchanger. 
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NOMENCLATURE 

English  Letter  Symbols 

A   =  Exchanger  heat-transfer  surface,  sq  m 

Am  =  Coefficient  of  the  mth  value  dimensionless 

Aq   =  Coefficient,  dimensionless 

Ai   =  1st  order  coefficient  to  be  multiplied  by  N-tu, 
dimensionless 

A2   =  2nd  order  coefficient  to  be  multiplied  by  N-tu2, 
dimensionless 

A3   =  3rd  order  coefficient  to  be  multiplied  by  N-^3, 
dimensionless 

A4   =  4th  order  coefficient  to  be  multiplied  by  N-tu4, 
dimensionless 

A5   =  5th  order  coefficient  to  be  multiplied  by  N-^5, 
dimensionless 

a    =  Exchanger  heat-transfer  surface,  sq  m/m 

C    =  Capacity  rate,  W/°K.   Also  designates  dimensionless 
arbitrary  constant 

Cpc  =  Specific  heat  at  constant  pressure  of  cold  fluid, 
J/kg°K 

Cph  =  Specific  heat  at  constant  pressure  of  hot  fluid, 
J/kg°K 

D   =  Empirical  value  of  effectiveness  (computer  generated) , 
dimensionless 

e    =  Error.   Also  used  as  the  exponential  function 

F    =  Logarithmic  mean  temperature  difference  correction 
factor,  dimensionless 

L   =  Exchanger  length,  m 

N   =  Number  of  effectiveness  empirical  data  points  used  to 
determine  a  curve  for  R,  dimensionless 


n   =  Number  of  tube  passes,  dimensionless.   Also  number  of 
equations,  dimensionless 

nc   =  A  related  Ntu  per  unit  length  cold  side,  m~l 

n^   =  A  related  Ntu  per  unit  length  hot  side,  m~l 

N-^u  =  Number  of  transfer  units,  dimensionless 

P   =  Temperature  group,  dimensionless 

q   =  Total  rate  of  heat  transfer,  W 

qmax=  Maximum  total  rate  of  heat  transfer,  W 

R   =  Capacity  rate  ratio,  dimensionless 

S   =  Temperature  group,  dimensionless 

Sr   =  Sum  of  the  squares  of  the  residuals,  dimensionless 

T   =  Hot  fluid  temperature,  °C 

Tpi  =  Particular  integral,  dimensionless 

Tj_   =  Hot  fluid  temperature  in,  °C 

T2   =  Hot  fluid  temperature  out,  °C 

ti   =  Cold  fluid  temperature  in,  °C 

t2   =  Cold  fluid  temperature  out,  °C 

ta   =  Cold  fluid  temperature  1st  pass,  °C 

tab  =  Cold  fluid  temperature  between  1st  and  2nd  passes 

tb   =  Cold  fluid  temperature  2nd  pass,  °C 

t^c  =  Cold  fluid  temperature  between  2nd  and  3rd  passes 

tc   =  Cold  fluid  temperature  3rd  pass,  °C 

tC(j  =  Cold  fluid  temperature  between  3rd  and  4th  passes 

t^   =  Cold  fluid  temperature  4th  pass,  °C 

tde  =  Cold  fluid  temperature  between  4th  and  5th  passes 

te   =  Cold  fluid  temperature  5th  pass,  °C 
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U   =  Overall  heat  transfer  coefficient,  W/m2  -  °C 

W   =  Mass  flow,  kg/sec. 

x   =  length  coordinate,  m.   Also  used  to  represent 
a  constant  value  in  a  sequence 

y   =  Sum  of  mth  degree  polynomial,  defined  by  eq.  (5.1), 
dimensionless 

Greek  Letter  Symbols 

a    =  Root  of  auxiliary  differential  equation,  1/m 

e    =  Exchanger  effectiveness,  dimensionless 

X    =  Combination  of  variables  defined  by  equation  (3.12), 
dimensionless 

£  =  Summation,  dimensionless 

a2  =  Variance,  dimensionless 

8m  =  Mean  temperature  difference  for  exchanger,  °C 

3  =  Indicates  partial  derivative,  dimensionless 

3  =  A  combination  of  terms  defined  by  eq.  (3.23) 

Subscripts 

c   =  Cold  fluid 

h   =  Hot  fluid 

i.j.k  -  Values  in  a  sequence 

m   =  Degree  or  order,  an  exponent 

1  =  inlet 

2  =  outlet 
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I.   INTRODUCTION 

In  many  engineering  applications,  it  is  desirable  to 
transfer  thermal  energy  between  fluids  at  different 
temperatures.   This  action  may  be  economically  accomplished 
with  the  aid  of  a  heat  exchanger. 

The  primary  function  of  a  non-direct  contact  heat 
exchanger  is  to  provide  two  paths  of  flow,  one  for  the  hot 
fluid  and  one  for  the  cold,  by  means  of  which  heat  can  be 
transferred  through  walls  separating  the  fluids. 

The  quantity  of  heat  transferred  is  governed  by  three 
factors:   1)  the  amount  and  nature  of  the  heat  transfer 
surface  to  which  the  two  fluids  are  exposed,  2)  the  overall 
heat  transfer  coefficient,  and  3)  the  mean  temperature 
difference  across  the  intervening  wall  from  one  fluid  to  the 
other. 

Numerous  flow  arrangements  and  geometric  designs  are 
possible.   One  of  the  most  common  types  of  heat  exchangers 
is  the  shell  and  tube  type.   Here  one  fluid  flows  through 
the  tubes.   The  tubes  are  enclosed  in  a  shell  with 
provisions  for  the  other  fluid  to  flow  through  the  spaces 
between  and  exterior  to  the  tubes.   Overall  flow  may  be 
either  parallel  or  counter,  but  in  most  cases,  the  latter  is 
preferred.   The  reason  for  this  will  be  seen  later.   One  of 
the  most  common  flow  arrangements  used  in  the  shell  and  tube 
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heat  exchangers  is  a  combination  of  the  two,  the  parallel- 
counterflow  exchanger.   Heat  exchangers  with  one  shell  pass 
and  n  tube  passes  are  therefore  referred  to  as  1-n  heat 

exchangers. 

One  example  of  this  flow  arrangement  is  the  one 
shell-pass  two  tube  pass  exchanger  as  shown  in  Figure  1.1. 


Baffle 


t 


I' 


t 


Figure  1.1   A  1-2  Heat  Exchanger  (One  Shell  Pass  and  Two 
Tube  Passes) 

This  type  of  exchanger  is  configured  so  that  all  of  the  tube 
side  fluid  flows  consecutively  through  half  of  the  total 
number  of  tubes:   The  baffles  in  the  shell  serve  to  "mix" 
the  shell  fluid  and  to  maintain  fluid  flow  at  right  angles 
to  the  tube  passes  in  order  to  obtain  a  maximum  "shell 
side"  heat  transfer  coefficient. 
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O'Hare  [Ref.  1]  has  indicated  that,  to  date,  much  work 
has  been  done  on  finding  the  logarithmic  mean  temperature 
difference  for  heat  exchangers  with  an  even  number  of  passes 
but  little  has  been  done  with  exchangers  having  an  odd 
number  of  passes.   From  this  basic  statement,  he  proceeded 
with  a  theoretical  examination  of  the  1-3  heat  exchanger 
which  resulted  in  some  startling  conclusions.   As  a  follow 
up  to  his  work,  this  thesis  will  parallel  and  follow  the 
format  of  O'Hare  [Ref.  1]  for  the  1-5  parallel  counter  flow 
exchanger.   One  arrangement  of  the  1-5  exchanger  is 
displayed  in  Figure  2. 
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Figure  1.2   A  1-5  Heat  Exchanger  (One  Shell  Pass  and  Five 


Tube  Passes) 
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II.   THE  DEVELOPMENT  OF  THE  EFFECTIVENESS  METHOD 

A.   HISTORICAL  DEVELOPMENT 

Nagle  [Ref.  2:   pp.  604-609],  in  1931,  credited  Davis 
[Ref.  3]  with  a  simplified  method  for  computing  actual  temp- 
erature differences  between  two  heat-interchanging  streams 
which  depart  from  true  counter  or  concurrent  (parallel) 
flow.   This  is  now  the  familiar  "F  factor"  method  which 
expresses  the  actual  mean  temperature  difference  8ra  in  q  = 
UA8m  as  a  fraction  F  of  the  counterflow  logarithmic  mean 
temperature  difference,  LTMD,  6rac  via  9m  =  F0mc. 

The  example  of  initial  interest  was  the  1-2  exchanger 
with  a  single  shell  pass  and  two  continuous  tube  passes  in 
counter  and  concurrent  flow  with  it.   The  method  involved 
derivation  of  the  actual  temperature  difference  for  the  flow 
pattern  and  formed  the  ratio  F  =  9m/9mc.   This  familiar  LMTD 
correction  factor  was  plotted  conveniently  as  functions  of 
the  effectiveness,  e,  and  the  capacity  rate  ratio  R  with  R 
as  a  parameter.   These  mean  temperature  difference 
correction  charts  are  available  for  many  flow  arrangements 
[Ref.  4:   pp.  829-833  and  Ref.  5].   The  effectiveness,  e 
(often  called  P  or  S),  is  always  the  cold  fluid 
effectiveness  and  R  is  always  the  capacity  rate  ratio  of 
cold  fluid  to  hot  fluid. 
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Nagle  detailed  assumptions  and  derivations  for  the  1-2, 
1-4  and  1-6  exchangers.   The  F  factors  were  obtained  by 
Nagle  through  graphical  integration  and  were  accompanied  by 
the  comment  that  F  factors  for  the  1-2  exchanger  could  be 
applied  with  negligible  error  to  1-4  and  1-6  exchangers. 
Underwood  [Ref.  6:  pp.  145-148]  rederived  the  equations  of 
Nagle  for  1-2  and  1-4  exchangers  to  eliminate  the  need  for 
obtaining  F  factors  by  graphical  integration. 

Bowman  [Ref.  7:  pp.  541-544]  pointed  out  that  for  a 
very  large  or  infinite  number  of  tube  passes,  the  F  factor 
approached,  as  a  limit,  its  value  in  crossflow  with  both 
fluids  completely  mixed.   It  was  further  stated  that  even  at 
the  limit,  the  F  factors  were  only  1  to  2  percent  lower  than 
those  for  the  1-2  exchanger.   A  previous  paper  by  Kraus  and 
Kern  [Ref.  8]  did  not  confirm  the  generalization  that  1-n 
exchangers  differed  only  negligibly  from  the  1-2  exchanger 
although  this  lack  of  confirmation  was  obtained  on  an 
e  =  f(R,  N-(-u)  basis.   Moreover,  the  Kraus-Kern  work  does  not 
confirm  the  generalizations  on  an  F  =  f(R,  N|-u,  e)  basis. 

From  the  standpoint  of  usefulness  and  good  accuracy,  it 
is  essential  that  F  factors,  if  they  are  to  be  used  in 
preference  to  e  =   f(R,  N-{-u,  flow  arrangement),  be  obtained 
with  precision.   Plots  of  F  =  f(R,  e,  =  P  or  S)  [Ref.  4: 
pp.  829-833  and  Ref.  5]  show  that  the  curves  for  particular 
values  of  R  approach  infinite  slope  as  F  decreases.   While 
this  can  be  partially  alleviated  by  restricting  R  <  1.0  (a 
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constraint  used  in  the  e  =  f(R,  N-tu,  flow  arrangement) 
approach),  it  is  seen  that  small  errors  in  the  interpolation 
for  R  or  e  =  P  or  S  can  result  in  large  fluctuations  in  the 
value  of  F. 

In  a  comprehensive  paper,  Bowman,  Mueller  and  Nagle 
[Ref.  9:  pp.  283-294]  presented  graphs  of  F  factors  for 
shells  with  one  through  six  shell  passes  and  numbers  of 
continuous  tube  passes  respectively  double  the  number  of 
shell  passes.   In  view  of  the  earlier  references  to  Nagle 
and  Bowman,  it  should  be  noted  that  F  factors  were  computed 
for  the  1-2  exchanger  in  [Ref.  9:  pp.  283-294]  using  the 
equations  of  Underwood  [Ref.  6:  pp.  145-158]. 

Ten  Broeck  [Ref.  10:  pp.  1041-1042]  prepared  a  graph 
of  the  dimensionless  groups  now  known  as  e ,  R  and  N-^u  for 
the  1-2  exchanger.  Such  a  graph  had  the  added  versatility 
of  simplifying  the  calculation  of  performance  in  a  given 
exchanger  when  operating  at  conditions  different  for  those 
for  which  it  was  designed.   Kays  and  London  [Ref.  11: 
pp.  63-74]  prepared  similar  graphs  and  tables  of 
e  =  f(R,  N-£U,  flow  arrangement)  for  the  1-2  exchanger  and 
for  several  cases  of  crossflow  and  periodic  flow. 

In  the  foregoing  discussion,  it  has  been  observed  that 
there  are  two  methods  for  the  design  and  analysis  of  a  heat 
exchanger.   These  are  the  so-called  "F  factor"  and  "e-N-tu" 
approaches  and  it  is  to  be  noted  that  the  Heat  Transfer 
Literature  contains  many  references  of  design  equations  for 
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either  of  these  methods  for  exchangers  that  have  one  shell 
pass  and  even  numbers  of  tube  passes. 

It  was  O'Hare  [Ref.  1]  who  pointed  out  that  there  may  be 
a  virtue  in  having  an  odd  number  of  tube  passes  particularly 
in  a  marine  application.   Such  an  arrangement  allows  the 
tube  side  fluid  to  enter  one  side  of  the  exchanger  and  leave 
from  the  opposite  side  (see  Figure  1.2). 

With  the  possible  use  for  an  odd  tube  pass  heat 
exchanger  established,  a  literature  search  was  conducted  to 
uncover  any  previously  solved  cases  of  1-n  heat  exchangers 
with  n,  the  number  of  odd  tube  passes.   O'Hare  [Ref.  1], 
Fischer  [Ref.  12],  and  Stevens,  et.  al.  [Ref.  13]  have 
dealt  with  this  type  of  exchanger.   Although  Stevens, 
et .  al.  deal  with  multiple  crossflow  passes,  all  provide  the 
foundation  for  the  work  reported  here  as  well  as  additional 
work  pertaining  to  the  1-3  exchanger.   O'Hare's  investiga- 
tion centered  around  the  development  of  a  solution  for  the 
1-3  exchanger  in  the  ef f ectiveness-Ntu  framework  while 
Fischer  derived  an  equation  using  the  mean  temperature 
difference  for  the  1-3  exchanger.   Moreover,  Fischer 
expressed  his  results  in  terms  of  the  correction  factor  F 
and  treats  only  the  case  of  the  1-3  exchanger.   This  thesis 
will  expand  on  the  Ntu  method  of  thermal  heat  exchanger 
design  by  covering  the  1-5  heat  exchanger.   The  work  will 
treat  both  the  three  counterflow  and  two  parallel  flow 
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and  the  two  counterflow  and  three  parallel  flow  (1-5: 3P) 
cases  (see  Figures  2.1  and  2.2). 
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Figure  2.1   1-5 :3C  Five  Tube  Passes  -  Three  in  Counterflow 
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Figure  2.2   1-5 : 3P  Five  Tube  Passes  -  Three  in  Parallel  Flow 
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III.   AN  ATTEMPT  AT  AN  ANALYTIC  SOLUTION 

A.   EFFECTIVENESS  AS  A  FUNCTION  OF  CAPACITY  RATE  RATIO  AND 
EXCHANGER  SIZE  (Ntu) 

This  section  deals  with  an  investigation  into  the 

effectiveness,  e ,  of  a  one  shell  pass-five  tube  pass  heat 

exchanger,  in  which  e  compares  the  actual  heat  transfer 

capability  to  the  therraodynamically  limited,  maximum  possible 

heat  transfer  capability.   This  exchanger  heat  transfer 

effectiveness  is  given  by 


C(T      -T        ")  CCt        -t        ) 

_     q       hv  hot,  in   hot,  out '  .   cv  cold, out    cold^n' 

e'a      CCT      -t       )  c        nr               -  t 

Mmax    minv  hot, in   cold, in'  minv  hot, in   cold, in) 


where  Cm^n  is  the  smaller  of  the  C^  and  Cc  magnitudes  (the 
capacity  rates).   Thus,  e  possesses  the  significance  of 
effectiveness  of  the  heat  exchanger  from  a  thermodynamic 
point  of  view,  with  the  magnitude  of  the  effectiveness 
completely  defining  the  heat  transfer  performance.   In 
general  we  express  e  =  f(N^u,  R,  and  flow  arrangement)  and 
when  the  flow  arrangement  is  specified  or  understood,  it  is 
said  that  e  =  f(Ntu>  R) .   [Ref.  11:   pp.  14-26]. 

The  number  of  heat  transfer  units  N^u  is  a  nondimensional 
expression  of  the  "heat  transfer  size"  of  the  exchanger. 
When  N-£U  is  small  the  exchanger  effectiveness  is  low,  and 
when  N-£U  is  large,  e  asymptotically  approaches  the  limit 
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imposed  by  the  flow  arrangement  and  thermodynamic 
conditions.   From  inspection  of  the  definition  of  Ntu 

A 


Ntu  =  T^   "  TT7- 
min    min 


lin   v 


it  is  clear  that  the  overall  conductance  and  transfer 
area  affect  the  costs  of  attaining  a  high  value  for  N-tu, 
ergo  a  high  e.   The  capacity  rate  ratio,  R,  as  defined  by 


C  . 
p    min 

R  "  c — 

max 


is  simply  the  ratio  of  mass  flow  rate  times  specific  heat 
capacity  for  the  two  streams.   These  can  be  considered  as 
flow  stream  thermal-capacity  rates,  i.e.,  the  energy  storage 
rate  in  the  stream  per  unit  of  temperature  change.   [Ref. 
11:  pp.  14-26] 

The  attempt  taken  in  this  thesis  to  develop  a  closed 
form  solution  has  used  the  basic  fundamentals  of  heat 
transfer  as  well  as  the  foregoing  definition.   A  closed  form 
solution  for  e  was  sought  for  both  1-5  exchangers  with  one 
having  three  out  of  five  tube  passes  in  parallel  flow  and 
the  other  having  three  out  of  five  tube  passes  in  counter- 
flow.   The  analytical  approach  taken,  and  demonstrated  in 
this  section,  is  for  three  out  of  five  tube  passes  in 
counterf low. 
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B.   ANALYTICAL  DEVELOPMENT 

The  derivation  for  the  effectiveness,  e,  of  the  1-5 
exchanger  as  a  function  of  the  capacity  rate  ratio,  R,  and 
number  of  transfer  units,  N-^,  depends  on  several 
assumptions. 

(1)  The  overall  coefficient  of  heat  transfer,  U,  does 
not  vary  within  the  exchanger. 

(2)  The  specific  heat  of  both  hot  side  and  cold  side 
fluids  does  not  vary. 

(3)  Each  fluid  is  thoroughly  mixed,  that  is,  the 
temperature  of  both  hot  and.  cold  side  fluids  is 
uniform  over  any  cross  section. 

(4)  Steady  flow  conditions  are  maintained. 

(5)  Heat  losses  to  or  from  the  environment  are 
negligible. 

(6)  No  change  of  phase  takes  place;  all  heat 
transferred  is  sensible  heat. 

(7)  There  is  equal  heat  transfer  surface  in  each  pass. 
The  configuration  is  shown  in  Figure  3. 1  where  the  five 

tube  passes  are  designated  with  subscripts  a,  b,  c,  d  and  e. 
The  temperature  of  the  hot  (shell  side)  fluid  is  indicated 
by  upper  case  letters.   For  the  cold  (tube  side)  fluid, 
lower  case  letters  are  used.   The  subscript  1  always  refers 
to  the  fluid  inlet  and  the  subscript  2  always  refers  to  the 
fluid  outlet. 

With  Wfo  and  Cp^  designating  mass  flow  (kg/sec)  and 
specific  heat  (Joules/kg~°K)  of  hot  fluid  entering  at  T^  and 
leaving  at  T2  we  define  a  capacity  rate  for  the  hot  side 

ch  =  whcph 
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Figure  3.1   Three  Counter  Two  Parallel  Pass  Configuration 
for  Development  of  the  Sought  After 
Effectiveness  Relationship 
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In  similar  fashion  for  the  cold  side  (with  Wc  and  Cpc) 
entering  at  t^  and  leaving  at  t2,  we  have 

We  then  obtain  an  energy  balance  for  the  entire 
exchanger 

Ch(Tl  "  T2>  "  cc(t2  "  ti)  (3.1) 

Over  the  right  hand  side  of  the  exchanger  (Figure  3.1) 

Ch(Ti  -  T)  =  Cc(t2  "  te  +  td  -  tc  +  tb  -  ta)  (3.2) 
and  a  differentiation  gives 

ChdT  =  Cc(dte  -  dtd  +  dtc  -  dtb  +  dta)  (3.3) 

Across  dx,  with  a  (ra2/m) ,  the  surface  per  running  meter 
of  length  of  pass  so  that  A  =  5aL  is  the  total  surface  in 
the  exchanger,  we  may  write  the  heat  transferred  to  the 
element  dx  in  each  cold  pass. 

Ccdta  =   Uadx(T  -  ta)  (3.4a) 

Ccdtb  =  -Uadx(T  -  tb)  (3.4b) 

Ccdtc  =   Uadx(T  -  tc)  (3.4c) 

Ccdtd  =  -Uadx(T  -  td)  (3.4d) 

Cddte  =   Uadx(T  -  te)  (3.4e) 
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Here  it  should  be  observed  that  due  cognizance  has  been 
taken  of  the  direction  of  the  flow  in  each  cold  fluid  pass 
with  respect  to  the  positive  sense  of  the  length  coordinate, 
x,  and  U  is  the  overall  heat  transfer  coefficient  (W/m^-°C). 
With  eqs.  (3.4)  in  eq.  (3.3) 

ChdT  =  Ua(5T  -  ta  -  tb  -  tc  -  td  -  te)dx 

or 

where 


Ua 
h 


is  a  sort  of  Ntu  per  unit  length  for  the  hot  side. 
Now  differentiate  eq.  (3.5) 


d2T      „dT   ^a   ^b   ^c   ^d   ^e  . 

T~2~    nhCbdx    dx     dx    dx     dx    dx   ; 
dx 


and  with  eqs.  (3.4)  substituted 


7^  •   5nh  Hi  -  nc"h(T  "  la  +  *b  "  *c  +  *d  "  V<3-6> 
dx 


where 


Ua 

nc  "  C- 
c 


where  again  the  resemblence  of  nc  to  Ntu  can  be  noted 
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From   eq.     (3.2)    we    obtain 


J±   (T.    -   T)    -   t0   =   t.     -   t      -    t      +    t ,    -    t  (3.7) 

Cvly2bacde  v/ 

c 


and  with  eq.  (3.7)  put  into  eq.  (3.6) 


XT"  5%S=  "Vh  IT  ^  (^  "  T)  -  t2] 
dx  c 


or 


dx  c 


(3.8) 


where 


Rc  ~  ^c/ch 

is  the  capacity  rate  ratio  for  the  cold  side 
Notice  that 


and 


Ch   Ua   Ua   Ch    ,Ua,2 

ncnh  c-  =  c--cr-c-=(c-)  =m 

c    c    h    c     c 


1  cv, 

R    =  -   h 

nh  "  R  C 

c  c 


a  capacity  rate  ratio  for  the  hot  side.   Then,  algebraic 
adjustment  provides 
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dx  h  h 


which  is  a  linear,  non-homogeneous,  second  order 
differential  equation  with  constant  coefficients  having  a 
complementary  function 

a- x      a9x 

T   =  C-e    +  C0e  (3.10) 

C     1  2  v      ' 

where  C^  and  C2  are  arbitrary  constants  and  where 


5nh  ♦  1  „,  2  a  /  :  v,1/2 

al'a2   =  -2-   *  5   [25nh      "    4m<      tt        > ' 

h 


5nh        %    ,„,         4m      }    -   Rh      ^2 


+ 


[25    -  ^    (—5 -)] 


But 


2      "    2       l ""  2    v      R, 

n,  h 

h 


2 

m      m    (Ua)2    .    ^h^  =    (Ch   2    =   R   2   =      1 

nT2"         (C    )2         (Ua)2  ^c"  h  R~7 

he  c 


so    that 


or 


5n.  n.  1    -   R        1/2 

»1'«2    -TT    *    2^    t25    "   4Rh      (~V^)] 


n. 
«1 ,a0    =  7^    (5    ±    A)  (3.11) 


l'u2         T~ 
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where 


X  =  [25  -  4Rh(l  -  Rh)]1/2  (3.12) 


Designate  the  particular  integral  as  Tp±   and  by  the 
method  of  undetermined  coefficients  let  Tp^  =  P  so  that  in 
eq.  (9) 


1  -  Ru        t9 
m(-*)  P  =  m(Hi-  Tx) 
h  h 


This  makes 


V  ■  p=  t  r"  Ti>  It-^tJ 

h  h 


so  that 


*2  -  Vi 

Pi  -   !  _  n  (3.13) 


"h 


The  general  solution  to  eq.  (3.9)  is  the  sum  of  eqs.  (3.10) 
and  (3.13) 

cux       a„x    t9  -  R.  T1 

T(x)  =  Cie    +  C2e  *      +  \    _\  (3.14) 

h 

where  the  arbitary  constants,  C^  and  C2  are  evaluated  from 
conditions  at  x  =  0  and  x  =  L.   At  x  =  0,  T(x  =0)  =  T2  and 
at  x  =  L,  T(x  =  L)  =  T^.   When  these  are  inserted,  in  turn, 
into  eq.  (3.14),  one  obtains  a  pair  of  linear  algebraic 
equations  in  the  unknowns  C^  and  C2 
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T0  =  C1  +  C0  +  T  . 
2     1     2     pi 


a1  L       ot„L 

T,  =  C.e    +  C0e    +  T  . 
11         2         pi 


where  Tp^  is  given  by  eq,  (3.13). 

It  is  only  a  matter  of  algebra  to  show  that 


and 


(T-  -  T  . )  -  (T0  -  T  . )e 
C-  =  -i ^ -4 Ei: (3.15a) 

e     -  e 


a  i  L 

(Tn  -  T  . )e     -  (T-  -  T  . N 

C   «  — 2 Eil -Li Eii  (3.15b) 

2  a.L     a„L 

e     -  e 


It  is  easy  to  see  from  eq.  (3.1)  that 


ch   (t2  -  y 

h   '  ^    (T1  "  T2) 


so  that 

t2    -    ti    +    Rh(Tl    "    T2) 
Use    of    this    in    eq.     (3.13)    shows    that 


T 
Pi  !    _ 


\ 
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or 

t   -  is :  y* 

xpi  ■   1  -  Rjj  (3.16) 

indicating  two  alternative  forms  for  Tp^  given  by  eqs. 
(3.13)  and  (3.16). 

Insertion  of  eqs.  (3.13)  and  (3.16)  in  eqs.  (3.15)  for 
C^  and  C2  will  yield  after  some  algebra 


,T1  -  V    ,T2  -  V   a2L 
<  1  -  ft  >  "  (  1  -  fth)  e 
C-  =   5_ , (3.17a) 

e     -  e 


and 


,T2  -  *lx   alL    ,T1  -  V 
(  1  -  R.  >  e    "  (  1  -  R  > 
C9  =   % 2_  (3.17b) 

e     -  e 


Equation  (3.14)  is  an  expression  for  the  hot  side 
temperature  at  any  location  in  the  exchanger  in  terras  of  the 
extreme  temperatures,  t^,  t2,  T]_  and  T2. 

Next  take  eq.  (3.5)  and  set  it  equal  to  the  derivative 
of  eq.  (3.14)  noting  that  C]_ ,  C2  and  Tp^  are  all  known 
constants. 

dT  alx        a2x 

Hx  =  nh(5T  "  *a  "  *b  "  *c  "  ^  "  te)  =  alCle    +  a2C2e 

(3. 18) 

At  x  =  0,  where  T  =  T2 ,  ta  =  ti,  tb  =  tc  =  tt>c  >  *d   =    le   =    tde 
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dT 

dl  =  nh(5T2  "  H   "  2tbc"  2tde  )  =  *i?±    +  a2C2      (3.19) 


and  if  we  subtract  eq.  (3.4a)  from  eq.  (3.4c)  we  obtain 


dt  -  dt      TT 

a     c     Ua  ,        , 

t  -  t   =  "  JT   dx  =  "ncdx 
a    c       c 


which  can  be  integrated  using  C3  as  the  constant  of 
integration 

-n  x 

t   -  t   =  C„e  c 
a     c     3 

At  x  =  0  where  ta  =  t^  and  tc  =  t^c 

tl    "  tbc  =  c3 

or 

tbc  =  tl  "  C3 
In  addition  at  x  =  L,  ta  =  tab  and  tc  =  tcd  so  that 


-n  L 

t  .  -  t  ,  =  C„e  C 
ab     cd     3 


or 


r   =  tab  "  tcd  =  ft    -  t   )enc 

C3 =TT>  (tab    cd;e 

c 

e 


This  gives  a  relationship  between  tab  and  tbc 


t    =t-rt    -  t   )e  (3.20) 

tbc    Zl         {   ab    cd;e  v     ' 
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where  Nc  =  ncL  can  be  considered  as  the  total  number  of 
transfer  units  for  the  cold  side. 

The  same  procedure  may  be  employed  to  find  t^e* 
Subtract  eq.  (3.4c)  from  eq.  (3.4e) 


dt   -  dt      TT 

c     e  _   Ua  ,   _   „  - 
_____ _.  dx  -  -  ncdx 

c    e       c 


and  by  integration 


-n  x 

t   -  t   =  C,e 
c    e    4 


At  x  =  0  where  tc  =  tfcc  and  te  =  t^, 
*bc  ~  *de  =  c4 


or 


^de  ~  ^bc  ~  c4 
At  x  =  L  where  tc  =  tC(j  and  te  =  t2 


"ncL 
ted  -  t2  =  c4e 


and  again  with  ncL  =  Nc 


and 


C4  =  (tGd  ~  t_)e  ° 


*de  -  *bc  "  (^d  "  *2)eNc  (3-21) 
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As  before  substract  eq.  (3.4b)  from  eq.  (3.4d) 


dt.  -  dt  .  TT 

b     d  Ua 

— r z —  =  TT—  dx  =  n  dx 

b    d  c 


and  by  integrating  find  that 


ncx 


at  x  =  0,  tb  =  tbc  and  td  =  tde  so  that 

tbc  ~  tde  =  ^5 

and  at  x  =  L,  tb  =  tab  and  td  =  tcd.   This  shows  that 

i 
tab  "  tcd  =  c5e 


ncL 


and 


Nc 
tcd  =  tab  "  (tbc  "  tde)e  (3.22) 


where  Nc  =  ncL. 

Return  to  eq.  (3.18)  and  look  at  the  conditions  at  x  =  L 
where  ta  =  tb  =  tab,  tc  =  td  =  tcd,  te  =  t2  and  T  =  Ti- 
These  conditions  in  eq.  (3.18)  give 

V5T1  -  2tab  "  2tcd  "  V  =  -ici^ai   +  -2C2eCt2 

(3.23) 
where  again  we  remember  that  C^  and  C2  are  known  constants. 

Now  with  four  equations,  eq.  (3.20),  eq.  (3.21), 
eq.  (3.22),  and  eq.  (3.23)  and  four  unknown  variables  we 
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can  solve  for  the  four  variables,  tab»  tbc ,  tC(j,  and  te . 
Cramer's  rule  will  be  used  to  solve  for  two  of  them;  t^c  and 

^de* 

First  arrange  eqs.  (3.20)  to  (3.23)  in  standard  form 

with  the  unknowns  appearing  in  columns. 


"Nc 
ab    be        cd 


=  t-e 


-N, 


(3.20) 


-  t 


be 


+  t  ,e  °  +  t  ,    =  t0e 
cd       de     2 


(3.21) 


tab  "  *bc     +  tcd    +  *de   =  ° 


-ab 


+  tcd 


=  6 


(3.22) 
(3.23) 


1        a-  L        a2L 
where  6  =  1/2(51^  -  tg  -  -=-    (a1C1e    +  a2C2e    ^ 

h 

If 


AX  =  B  = 


-N, 


1   e 


0    -1 


1   e 


N, 


-1 

1 

0 

eN<= 

1 

-1 

eN= 

1 

0 

ab 


be 


cd 
ide. 


t-  e 


-N, 


t2e 


0 


N, 


The  solution  for  t^c  and  t^e    is 


det  A 


be 


det  A 


de 


'be     det  A 


de     det  A 


36 


where,  AlDC  and  A^e  are  the  matrices  obtained  by  replacing 
the  entries  in  the  second  and  fourth  columns  of  A  by  the 
entries  in  matrix  B. 

Compute  the  determinant,  det  A  of  matrix  A  (the  matrix 
of  coefficients) 


det  A   = 


-N, 


1   e 


0    -1 


1    e 


N, 


-1     0 


Nc 
e  c    1 


-1 


N, 


0 


det  A   =  -(e  c  +  2e   °) 


Since  det  A   *  0 ,  the  system  should  have  a  unique  solu- 
tion.  Continue  by  computing  det  k^c      and  det  A^g.   Thus 


det  Abc   = 


~Nc 
1    txe      -1     0 


0     t2e  c      e  c    1 


-1 


N, 


2Nr  Nc      -N~ 

=  e  c(3  -  2t2)  -  tx(e    +  2e   c) 
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and 


det  Ade   = 


-N, 


0  -1 

1  -e  c 
1      0 


-1 


-N, 


-1     t1e 


Nc        N, 
e      t2e 


0 


Therefore 


=  (e   c  +  1)(3  -  2t9)  -  t, (e  c  +  2e   c) 


det  A 


be 


be    det  A 


e   C(B  -  2t2)  -  t1(e  c  +  2e   c) 

N       =N 
-(e  c  +  2e   c) 


or 


and 


be     1 


Nc 
e  c(3  -  2t2) 


-2N 


(3.24) 


1  +  2e 


2N  N       -N 

det  Ade     (eJ  c  +  1)(3  -  2tg)  -  t1(e  c  +  2e   c) 


de    det  A 


N       -N 
-(e  c  +  2e   c) 


or 


de     1    v 


2N 
(e   c  +  1)(B  -  2t2) 

N       =N 
e  c  +  2e   c 


(3.25) 
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Then  with  eq.  (3.24)  and  (3.25)  in  eq.  (3.19)  and  some 
algebraic  manipulation 

(3  -  2t9)(2e  c  +  e   c) 

«1C1  +  «2C2  =  V5<T2  "  V  +  2< =2N )] 

(1  +  2e   c) 

(3.26) 

Equation  (3.26)  is  similar  to  the  one  derived  by  O'Hare 

[Ref.  1:  eq.  (22)].   At  this  point  he  continued  to  attempt  a 

closed  form  solution  for  e  as  a  function  of  R  and  N^u.   It 

is  to  be  noted,  however,  that  he  was  unable  to  obtain  the 

solution  because  of  the  presence  of  terms  deriving  from 

dimensional  parameters  such  as  n^  with  dimensions  m~l 

[Ref.  1:  pp.  46].  Since  n^  is  present  in  eq.  (3.26),  it 

must  be  concluded  that  a  further  attempt  at  an  analytic 

solution  for  the  effectiveness  of  the  1-5  exchanger  will 

suffer  a  similar  fate  and  be  also  unachieveable. 

This  conclusion  should  be  obvious  because,  when  dealing 

with  an  equation  derived  from  n  equations  with  n+1  unknowns, 

one  cannot  create  an  additional  equation  by  merely 

multiplying  one  of  the  n  equations  by  a  constant.   This 

leads  to  a  linearly  dependent  set  which  has  no  unique 

solution. 

Thus,  after  observing  that  O'Hare's  attempt  failed,  this 

thesis,  like  O'Hare's,  turns  to  a  numerical  analysis  to 

achieve  the  requisite  objective. 
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IV.   NUMERICAL  AND  COMPUTER  ANALYSIS 

As  O'Hare  discovered  with  the  1-3  heat  exchanger 
[Ref.  1],  an  analytic  solution  for  effectiveness  as  a 
function  of  the  capacity  rate  ratio,  R,  and  N-£U  was  also 
unattainable  for  the  1-5  exchanger.   Hence,  an  alternate 
approach  had  to  be  taken  and  this  involved  utilizing  a 
finite-difference  thermal  analyzer  combined  with  a 
linearizing  scheme  for  the  solution  of  a  set  of  nonlinear 
algebraic  equations  [Ref.  14].   This  thermal  analyzer  has 
been  employed  to  provide  solutions  for  temperature  that  lead 
to  effectiveness,  e,  value  for  the  1-5 :3C  and  1-5: 3P  heat 
exchangers  as  a  function  of  R  and  N-(-u. 

A.   THERMAL  ANALYZER  TVSSI 

The  computer  program  TVSSI  used  in  the  development  of 
the  e  =  f(R,  N-tu)  graphs  for  the  1-5  heat  exchanger  [Ref.  1: 
Appendix  A],  was  adapted  from  TVSS2  listed  by  Kern  and  Kraus 
[Ref.  14:   Appendix  C].   This  adaptation  consisted  of  modi- 
fying TVSS2  to  compute  solutions  in  the  SI  system  and  accept 
an  input  file  specifically  designed  for  the  1-5  heat  exchan- 
ger.  This  program  also  had  to  be  adapted  to  the  IBM  3033  AP 
system.   The  program  utilizes  the  Cholesky  decomposition 
scheme,  and  because  of  the  linearization  of  the  radiation 
terms  (not  used  in  this  study)  the  program  is  iterative. 
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B.   DEVELOPMENT  OF  THE  1-5 :3C  and  1-5 :3P  HEAT-EXCHANGER 
MODELS 

As  indicated  in  the  foregoing,  TVSSI  uses  an  input  file 
to  generate  a  solution  for  temperatures  that  will  lead  to 
the  sought  after  effectiveness  as  a  function  of  R  and  N-(-u. 
This  input  file  is  created  by  a  program  which  requires  an 
input  of  given  capacity  rates  (hot  and  cold) ,  overall 
coefficient,  surface  area  and  inlet  stream  temperatures  (hot 
and  cold)  values. 

This  required  the  development  of  a  program  which  created 
such  a  file.   O'Hare  [Ref.  1:  pp.  53-55]  first  developed  a 
program  using  the  1-4  heat  exchanger  that  would  create  an 
input  file  for  TVSSI.   He  compared  the  effectiveness  from 
the  computer  generated  results  and  a  known  analytic  solution 
of  effectiveness  developed  by  Kraus  and  Kern  [Ref.  8].   Once 
O'Hare  established  confidence  in  his  results,  he  proceeded 
to  develop  a  program  which  would  create  an  input  file  for 
the  1-3  exchanger.   In  this  study,  a  similar  technique  was 
used  to  develop  a  computer  program  for  the  1-5 :3C  and  1-5 :3P 
heat  exchangers.   These  programs  are  called  NTU53C  and 
NTU53P  (see  Appendices  A  and  B)  and,  in  them,  the  following 
parameters  and  techniques  were  utilized. 

1.  300  nodes  were  used. 

2.  The  initial  temperature  for  the  computer  to  begin  the 
iterative  process  was  set  at  200°C. 

3.  An  eventual  accuracy  of  .05  between  the  final  and  next 
to  last  iterations  was  used. 
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4.  A  radiation  coefficient  convergence  factor  of  0.66667 
between  iterations  was  used.   This  is  required  by  the 
program  even  if  radiation  is  not  considered  as  a 
node. 

5.  The  maximum  number  of  iterations  that  the  computer 
was  allowed  to  perform  was  set  at  12. 

6.  A  damping  factor  of  .8  was  set  as  an  initial  damping 
based  on  the  number  of  non-linear  terms  in  all  of  the 
node  equations. 

With  the  values  of  C^,  Cc ,  U,  A,  T^  and  ti  set,  the 
input  is  set  and  an  input  file  for  TVSSI  was  generated.   In 
this  file  all  node  equations  and  internode  conductance 
values  were  determined.  • This  program  specified  the  nodes 
that  interact  with  each  other  and  the  methods  by  which  the 
thermal  interaction  takes  place  such  as  in  the  forced 
convection  and  fluid  flow  modes. 

Both  programs  make  use  of  the  fact  that  each  term 
in  a  node  equation  shows  three  things.   The  first  is  the 
node  that  is  coupled  for  heat  flow  with  the  node  in 
question.   The  second  is  the  method  of  heat  flow  between  the 
nodes.   In  this  case, forced  convection  and  fluid  flow  are 
used.   Finally,  the  node  equation  shows  the  magnitude  of  the 
internode  heat  flow.  Here  all  the  pieces  of  information  are 
collected  and  presented  for  use  by  TVSSI  as  an  input  file 
with  all  items  in  the  proper  format. 

It  is  the  output  values  of  both  NTU53C  and  NTU53P  (a 
typical  result  is  shown  in  Appendix  C)  that  are  used  by 
the  thermal  analyzer  (as  an  input  file)  to  determine 
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the  temperatures  T2  and  t2  for  the  specific  set  of  given 
initial  parameters  C^,  Cc,  R,  A,  T^,  and  t]_. 

C.   SCOPE  OF  COMPUTER  ANALYSIS 

With  the  development  of  NTU53C  and  NTU53P  computer 
programs,  the  files  that  they  generated  could  then  be  used 
to  solve  for  temperatures  which  in  turn  could  yield  effec- 
tiveness values  for  a  given  R  and  N-(-u.   Of  course,  it  must 
be  realized  that  many  computer  runs  are  required  to  generate 
enough  data  to  ensure  confidence  in  the  results  which  cover 
a  wide  range  of  capacity  rate  ratios  and  N^u.   values. 

To  efficiently  expedite  the  computer  task,  the  Multiple 
Virtual  System  (MVS )  with  Job  Entry  Subsystem  and  Networking 
(JES3)  was  employed.   The  MVS  coupled  with  JES3  is  more 
commonly  referred  to  as  batch  processing.   Based  on  trial 
and  error,  it  was  determined  that  in  order  to  build  a  solid 
data  base,  eleven  different  values  for  effectiveness  were 
needed  to  best  represent  a  particular  value  of  R.   Thus,  a 
range  of  R  from  R  =  0.01  to  1.0  was  used.   Due  to  the 
enormous  amount  of  computer  time  required  to  solve  a  300 
node  problem,  R  increments  of  0.05  were  used  for  this  study 
and  this  thesis  reports  on  this  basis.  In  all,  42  curves  for 
both  the  1-5 :3C  and  1-5 :3P  heat  exchangers  were  required. 
With  this  change  in  the  increments  between  R  values 
(0 'Hare's  work  showed  that  an  increment  in  R  of  0.01  was  not 
necessary),  a  linear  interpolation  was  found  to  be  accepta- 
ble to  less  than  0.02  percent  error. 
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TVSSI  was  slightly  modified  in  accordance  with  the 
appropriate  guidelines  of  the  job  control  language  (JCL) 
needed  to  run  on  the  batch  processing  system.   These 
modifications  were  few  and  were  needed  only  at  the  beginning 
and  end  of  TVSSI.   The  modified  version  of  TVSSI  has  been 
called  TVCOUNT  with  changes  shown  in  [Ref.  1:  Appendix  F]. 
It  was  TVCOUNT  that  was  then  used  to  activate  TVSSI. 

It  also  became  necessary  to  modify  the  two  input  file 
programs  NTU53P  and  NTU53C  so  that  they  needed  to  be 
compiled  only  once.   They  were  then  loaded  in  a  library 
file  to  be  used  when  called  by  another  program.   New 
programs  utilizing  the  batch  system  were  written  that  could 
easily  be  loaded  with  the  appropriate  input  data  for  a 
specific  R  value.   These,  which  are  referred  to  as  "sister 
programs,"  were  used  to  go  from  the  library  file  to  TVSSI 
and  cause  TVSSI  to  be  executed  eleven  times  under  TVCOUNT 
covering  the  desired  range  of  N-^  for  a  specific  R  value. 
The  modified  input  files  and  their  associated  "sister 
execution  programs  are  found  in  Appendices  D  through  G. 

A  flow  chart  of  how  all  of  the  foregoing  is  accomplished 
is  provided  by  O'Hare  [Ref.  1:   Figures  4.4  and  4.5].   It  is 
noted  that  in  these  figures,  TVSSI  is  referred  to  as  TVSSIA 
through  TVSSIV.   These  are  the  same  programs  as  TVSSI  but 
for  bookkeeping  purposes  by  the  computer  they  are  labeled  A 
through  V. 
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D.   INTERPRETATION  OF  GRAPHS 

Once  all  the  data  had  been  collected  and  the 
effectiveness  for  a  given  R  and  N-^  values  had  been 
computed,  plots  of  effectiveness  as  a  function  of  N-^u  for 
the  entire  range  of  R  values  were  plotted  for  both  parallel 
and  counterflow  heat  exchangers.   These  plots  are  shown  in 
Appendices  H  and  I. 

It  may  be  observed  that  for  a  particular  value  of  R,  the 
plot  of  effectiveness  as  a  function  of  Ntu  indicates  that 
the  1:5-3C  exchanger  outperforms  the  1-4,  1-3: 2P  and  the 
1-5: 3P  exchangers  but  not  the  1-3: 2C  exchanger.   This  may  be 
explained  by  the  quite  apparent  fact  that  the  1-5: 3C 
exchanger  has  one  more  counterflow  pass  than  the  1-4 
exchanger.   Indeed,  this  argument  can  be  extended  to  the 
fact  that  the  1-5: 3C  has  one  more  parallel  flow  pass  than 
the  1-3: 2C.   Figure  4.1.  displays  the  graphical  picture  to 
support  this  contention. 


45 


■ 

o 

II 

CO 

DC 
O 


CO 

> 

CO 
CO 

LU 

z 

LLI 
> 

O 
4J 

J- 
4- 
LU 


o     o      a.  a. 

H       •»         «•    « 


\\lll 

• 

\\\r, 

1 

■ 

\ 

■ 

\\ 

V 

• 

\ 

■ 

\  x 

> 

\ 

■ 

■ 

* 

3 
z 


SS3N3ALLD3JJ3 


Figure  4 . 1 


Comparison  of  Analytical  (1-4)  and  Computer  Results 
(1-3-.2C  and  1-3:2P)  to  Computer  Results  (1-5:3C 
and  1-5:3P)  at  R  =  0.5 
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V.   POLYNOMIAL  REGRESSION 

A.   DEVELOPMENT  OF  POLYNOMIAL  EQUATIONS 

The  empirical  data  obtained  for  the  two  1-5  heat 
exchangers  was  designed  to  cover  an  extensive  range  of  R 
values  varying  from  0.01  to  1.0  in  increments  of  0.05.   As 
discussed  earlier  in  Section  IV,  the  data  obtained  at  a 
specific  value  of  R  is  the  computer  evaluated  result  of  the 
effectiveness,  for  an  associated  N-j-u  value.   With  this 
accomplished,  it  then  becomes  possible  to  graph  separate 
curves  for  each  of  the  different  R  values  as  shown  in 
Appendices  H  and  I.  Through  a  polynomial  regression 
technique,  as  discussed  in  this  section,  it  is  also  possible 
to  develop  implicit  equations  for  the  curves  with  e  =  f(N^Ul 
R).   It  is  also  apparent  from  an  inspection  of  the  graphical 
representation  of  the  empirical  data  in  Appendices  H  and  I, 
that  the  curves  conform  to  a  high  degree  polynomial. 
However,  further  analytical  investigation  is  needed  to 
ascertain  the  exact  order  of  the  polynomial  terms.   This 
investigation  will  not  only  lead  to  the  order  of  the 
polynomial,  but  to  the  specific  equation  for  each  curve. 

By  use  of  polynomial  regression,  the  least-squares 
method  can  be  readily  extended  to  best  fit  the  data  to  the 
mth-degree  for  the  polynomial 

y  =  A0  +  A]_x  +  A2x2  +  ...  Amxm  (5.1) 
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with  the  error  defined  by 

ei  =  Dj_  -  yi  =  Dj_  -  A0  -  Aj_x  -  A2Xi2  -  ...  -  AmXim 

where  D^  represents  the  empirical  data  value  corresponding 
to  x^  ,  Xi  being  free  of  error. 

The  objective  is  to  minimize  the  sum  of  the  squares  of 
the  residuals,  Sr, 


m     2    m  2  m  2 

S   =  Z      e.   =  Z    (D.  -  A_  -  A,x,  -  A0x0   +  ...  A  x  )  (5.2) 
r    .  _1   l        i    0     11    22         m 


Because  at  a  minimum,  the  partial  derivatives  3Sr/3AQ, 
3Sr/3Ai  ...  3Sr/3Am  vanish,  after  taking  the  derivative  of 
Sr  with  respect  to  each  of  the  coefficients  of  the  poly- 
nomial, it  can  be  seen  that 


=  0  =  -2  E  (D.  -  A_  -  A.x.  -  A_x.   -  ...  -  A  x.m) 


3  A 


as  9 

=  0  =  -2  Z    x.(D.  -  An  -  A,x.  -  A0x.   -  ...  -  A  x.m) 


3A     ~     ~    "iv  i    "o     1  i     2  i  mi 


=  0  =  -2  Z    x.  (D.  -  A_  -  A.x.  -  A0x.   -  ...  -  A  x.m) 


3A     ~         «±  v^i    ..Q    „1„i    „2„i  ^^ 


3S 

-^  =  0  =  -2  Z    x.m(D.  -  An  -  A,x.  -  A_x.   -  ...  -  A  x.m) 
3A  l  v  l     0     1  l     2  l  mi 

m 
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Then  by  dividing  by  -2  and  rearranging  we  obtain 


A0N  +  Ai  E  Xi  +  Axxi2  +  ...  +  Am  E  Xim  =  I    D± 

A0  E  x±    +  Ai  E  x±2    +  A2  E  Xi3  +  ...  +  Am  E  x^m+l  =  z  x^ 

A0  E  xt2  +  Ax  E  Xj^3  +  A2  E  Xi4  +  ...  +  Am  E  Xim+2  =  E  xi2Di 


A0  E  Xim  +  Ai  E  Xim+i  +  A2  E  x-l 


m+1 


m+2 


+  . 


+  Am  E  x±2m   =    E  xj^Di 


where  all  summations  are  from  i=l  through  n.   All  of  the 
foregoing  m+1  equations  are  linear  and  have  m+1  unknowns: 
Aq,  Ai ,  A2 ,  ...  Am.   The  coefficients  of  the  unknowns  can  be 
calculated  directly  from  the  observed  data.   Thus,  the  pro- 
blem of  determining  a  least-squares  polynomial  of  degree  m 
is  equivalent  to  solving  a  system  of  m+1  simultaneous  linear 
equations.   Putting  the  equations  in  matrix  form  yields 


E  x 


E  x 


E  xt2 


E  Xi2    E  xt3 


E  xt2 
E  xt3 
E  x-i4 


E  x- 


E  x 


E  x 


E  Xjm    E  Xim+1   E  Xim+2   E  Xim+3  ...Ex 


E  x-im 


E  Xim+1 


E  x 


m+2 


.2m 
■l 


A   = 


E  D- 


E  x^ 
E  xi2D± 


E  ximDi 


[Ref.  15:  pp.  302-309  and  Ref.  16:  468-474]. 

49 


From  this  point  on,  one  finds  that  it  is  best  to  use  a 
computer  to  assist  in  solving  the  simultaneous  equations  and 
this  will  also  help  alleviate  any  ill-conditioning  that  may 
otherwise  occur.   An  existing  "curvefit"  program  available 
through  NON-IMSL  [Ref.  16]  was  used  although  some  modifica- 
tions were  made  to  the  original  program  to  best  accommodate 
the  goals  of  this  work. 

To  determine  the  order  of  polynomial  that  should 
eventually  be  used,  one  increases  the  degree  of  the 
approximating  polynomial  as  long  as  there  is  a  statistically 
significant  decrease  in  the  variance  a*,    which  is  computed 
by 

2 

2     Z  ei 

a      =  ct r  (5.3) 

N  -  m  -  1  v    J 

In  otherwords,  the  selection  of  the  optimum  degree  poly- 
nomial is  contingent  upon  a  decreasing  variance  and  once  the 
variance  begins  to  increase,  the  degree  of  the  polynomial 
becomes  too  high.   For  all  cases,  it  was  found  that  the 
e  -  N-j-u  developed  curves  are  of  the  5th  order. 

As  shown  in  Figures  5.1  and  5.2  the  computed  values  of 
effectiveness  vs.  Ntu  for  R  =  0.1,  0.5  and  1.0  for  both 
flow  arrangements,  (1-5: 3P)  and  (1-5: 3C),  have  been  graphed 
and  fitted  by  a  5th  degree  polynomial.   Because  all  computed 
values  for  effectiveness  follow  a  predictable  trend, 
only  a  sample  of  the  data  covering  the  whole  range  of 
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Figure  5.1   1-5 :3C  Data  Fit  by  a  5th  Order  Polynomial 
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Figure  5.2   1-5:3P  Data  Fit  by  a  5th  Order  Polynomial 
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values  of  R  have  been  shown.   It  is  clear  that  the  graphical 
interpretation  strongly  backs  what  is  known  numerically 
from  the  polynomial  regression  technique.   Thus,  the 
relationship  for  e  =  f(Ntu,  R)  can  be  found  explicitly  from 


e  = 


A5Ntu5  +  A4Ntu4  +  A3Ntu3  +  A2Ntu2  +  AlNtu  +  A0  (5.4) 


where  the  coefficients  A5,  A4,  A3,  A2,  A^,  and  Aq  relating 
to  a  specific  value  of  R  are  found  in  Tables  1  and  2  for  the 
(1-5:3C)  and  (1-5:3P)  configurations. 

In  order  to  better  utilize  the  coefficients  found  in 
Tables  1  and  2,  an  interactive  program  has  been  developed 
using  equation  (5.4)  and  the  1-5  heat  exchanger  polynomial 
approximation  coefficients.   This  program  was  written  for 
the  IBM  3033  using  the  FORTRAN  VS  compiler  (See  Appendix  J). 
By  entering  the  values  for  NTU,  R  and  the  type  of  heat 
exchanger  the  effectiveness  values  are  readily  obtained 
directly  or  by  linear  interpolation. 

An  example  of  how  to  work  with  the  effectiveness  by 
using  the  Equation  (5.4)  follows. 

B.   NUMERICAL  EXAMPLE 

1.   Statement  of  Problem 

A  1-5: 3C  heat  exchanger  uses  sea  water  (Cp  =  4.045 
kJ/kg-°C)  at  an  initial  temperature  of  20°C  to  cool  engine 
oil,  flowing  at  200  kg/s  from  60°C  to  40°C.  The  sea  water 
flow  rate  is  80  kg/s.   If  the  overall  heat  transfer 
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coefficient  is  40  W/m2oC,  how  much  surface  area  is  required 
in  the  exchanger? 
2.   Solution 

At  an  average  oil  temperature  of  50°C,  the  specific 
heat  is  Cp  =  2.004  kJ/kg-°C  and  an  energy  balance  gives 

(WCpAT)sea  water  =  (^pA^)engine  oil 
(80  kg/s)(4.045  kJ/kg-°C)(t2  -  20)°C 

=  (200  kg/s)(2.004  kJ/kg-°C)(60  -  40)°C 


200(2.004)(20)  =  44  gOC 
t2    ZiJ    +    10(4.045)      44,8  L 


The  number  of  transfer  units  is  given  by  N-£U  = 
UA/Cm^n.   To  determine  Crain  we  must  compare  the  products  of 
mass-flow  rate  and  specific  heat  for  sea  water  and  oil. 

(wcp)sw   "  (80  kg/s)(4.045  kJ/kg-°C) 

=  323.6  kW/°C 

(wCp)oil  =  (20°  kg/s)(2.004  kJ/kg-'C) 

=  400.8  kW/°C 

Therefore  Cmin  =  323.6  kW/°C. 

The  effectiveness  is  a  function  of  R  and  N-tu.   The 
capacity  rate  ratio,  R,  is 
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R   = 


C  . 

mm 

C 
max 


323.6 
400.8 


=   0.8074 


Because  Cc  =  Cm^n,  the  required  effectiveness  is 


*2  -  h 
Tl  -  rl 


44.8  -  20 
60  -  20 


=   0.6193 


We  must  now  go  to  Table  1  for  the  1-5 :3C  arrangement  with 
R  =  0.8  and  R  =  0.85  and  find  the  coefficients 


A0 
Al 
A2 

A3 

A4 

A5 


R  =  0.8 

3.63419  x  10-3 

0.93822 
-0.67053 

0.27531 
-0.60558  x  10_1 

0.54393  x  10-2 


A0 
Al 
A2 

A3 

A4 

A5 


R  =  0.85 

3.82659  x  10"3 

0.93407 
-0.67841 

0.28109 
-0.62088  x  10"1 

0.55850  x  10-2 


A  trial  and  error  solution  (assuming  values  of  Ntu  and  by 
interpolating)  gives  a  result 


N 


tu 


=   2.352 


for  which 


A  =  Ntu  Crain/U 


A  =  19.03  m^ 
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3.   Observations 

The  primary  observation  made  here  is  that  by  using 
the  5th  order  polynomial  equation  with  the  appropriate 
coefficients  found  in  Table  1  or  2,  an  accurate  value  for 
effectiveness  can  be  found  thus  allowing  one  to  solve  for 
a  wide  variety  of  heat  exchanger  parameters. 

By  comparing  the  value  for  effectiveness  just 
computed  with  that  of  the  1-5: 3P,  1-3 :2P  [Ref.  1]  and  1-4  to 
1-12  even  pass  exchangers  [Ref.  8],  the  1-5: 3C  has  a  higher 
ef f-ectiveness  than  all  (see  Table  3). 

From  this  observation  it  is  evident  that  the  1-5 :3C 
exchanger  outperforms  not  only  the  1-4  to  1-12  even  pass 
arrangements  but  also  its  counterpart,  the  1-5: 3P  exchanger 
for  these  given  values  by  at  least  0.614%. 
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TABLE  3 
COMPARISON  OF  EXAMPLE  RESULTS 
FOR  Ntu  =  2.352  AND  R  =  0.8074 


NUMBER  OF 
TUBE  PASSES 

EFFECTIVENESS 

3:2P 

.5990 

4 

.6155 

5:3C 

.6193 

5:3P 

.6067 

6 

.6145 

8 

.6141 

10 

.6140 

12 

.6139 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

The  objective  of  this  thesis  was  to  provide  effective- 
ness values  for  the  1-5: 3C  and  1-5: 3P  heat  exchangers. 
Although  these  objectives  were  analytically  unachievable, 
they  became  attainable  through  a  fifth  order  polynomial 
approximation  from  numerical  data.   The  data  obtained  covers 
a  wide  range  of  capacity  rate  ratios  for  values  of  N-£U  from 
0  to  3.25.   With  the  knowledge  of  the  1-5  heat  exchanger 
type  (1-5: 3C  or  1-5: 3P),  size  (surface  area),  the  overall 
heat  transfer  coefficient  and  the  fluid  flow  rates,  the 
N-(-u  and  capacity  rate  ratios  may  be  computed  and  the 
exchanger  effectiveness  may  then  be  determined  from  the 
appropriate  coefficient  tables  or  the  effectiveness  program 
provided. 

It  is  indeed  obvious  that  the  effectiveness  parameter  is 
most  useful  in  performance  estimation  and  analysis  for  all 
types  of  heat  exchangers.   The  work  done  here  pertaining  to 
the  1-5  exchanger  is  very  helpful  in  constructing  a 
comprehensive  picture  of  the  performance  attainable  for  the 
1-2  consecutively  through  the  1-6  for  a  diverse  range  of 
known  conditions. 
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B.   RECOMMENDATIONS 

With  proof  in  this  thesis,  and  that  of  O'Hare  [Ref.  1], 
that  an  odd  number  of  tube  passes  per  shell  has  a  slightly 
better  effectiveness  when  more  than  half  the  tube  fluid 
flows  in  counterflow  to  the  shell  fluid,  it  is  recommended 
that  a  feasibility  study  be  conducted  with  regard  to  the 
structural  and  thermal  problems  that  may  arise  in 
manufacturing  and  designing  such  heat  exchangers. 
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APPENDIX  H 
1-5 :3P  EFFECTIVENESS  VS.  N    GRAPHS  AT  VARIOUS  R  VALUES 
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APPENDIX  I 
1-5 :3C  EFFECTIVENESS  VS.  N    GRAPHS  AT  VARIOUS  R  VALUES 


51 

■ 

CC 


CO 

> 

CO 
CO 
LU 


3 
O 

o 

z 


CO 
IU 

uj  S 
>£ 

to 
u. 

o 

O 

CO 


5 

UJ 
Li. 
LL 
LU 


-»0oc>«jo>dooo«5  snocoaoaa 


0000000000  o  o  oooooooo^ 

1  1    1    1   1    1    I    1    1   ■  1  1   1  1  1  I  I  1  1  1  1 

ccccccocccoccc  ccccccccecucccsc 

-1 

1 

\\ 

1  1  1  1 

1 

ill 

11 

\\\l 

III 

\\\\\ 

1 

w 

11 

I 

It 

it 

SI 

v\\w\\ 

■ 

\ 

11 

k 

^ 

%K 

^t® 

||l 

^~ 

' 

d 


CD 

d 


IN 

d 


© 
d 


SS3N3ALL03Jd3 


Figure  1.1 


1-5 :3C  Effectiveness  vs. 
0.01  to  1.0 


N    over  Range  of  R  from 


102 
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